Non-perturbative conserving approximations and Luttinger's sum rule 
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Weak-coupling conserving approximations can be constructed by truncations of the Luttinger- 
Ward functional and are well known as thermodynamically consistent approaches which respect 
macroscopic conservation laws as well as certain sum rules at zero temperature. These properties 
can also be shown for variational approximations that are generated within the framework of the self- 
energy-functional theory without a truncation of the diagram series. Luttinger's sum rule represents 
an exception. We analyze the conditions under which the sum rule holds within a non-perturbative 
conserving approximation. Numerical examples are given for a simple but non-trivial dynamical two- 
site approximation. The validity of the sum rule for finite Hubbard clusters and the consequences 
for cluster extensions of the dynamical mean-field theory are discussed. 

PACS numbers: 71.10.-w, 71.10.Fd 
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I. INTRODUCTION 

Continuous symmetries of a Hamiltonian imply the 
existence of conserved quantities: The conservation of 
total energy, momentum, angular momentum, spin and 
particle number is enforced by a not explicitly time- 
dependent Hamiltonian which is spatially homogeneous 
and isotropic and invariant under global SU(2) and U(l) 
gauge transformations. For the treatment of a macro- 
scopically large quantum system of interacting fermions, 
approximations are inevitable in general. Approxima- 
tions, however, may artificially break symmetries and 
thus lead to unphysical violations of conservations laws. 

Baym and Kadanoffi^ have analyzed under which cir- 
cumstances an approximation for time-dependent corre- 
lation functions, and for one- and two-particle Green's 
functions in particular, respect the mentioned macro- 
scopic conservation laws. They were able to give cor- 
responding rules for a proper construction of approxima- 
tions, namely criteria for selecting suitable classes of dia- 
grams, within diagrammatic weak-coupling perturbation 
theory. Weak-coupling approximations following these 
rules and thus respecting conservation laws are called 
"conserving" . Frequently cited examples for conserving 
approximations are the Hartree-Fock or the fluctuation- 
exchange approximation j 1 ' 3 ' 4 

Baym2 has condensed the method of constructing con- 
serving approximations into a compact form: A con- 
serving approximation for the one-particle Green's func- 
tion G is obtained by using Dyson's equation G = 
1/(Gq 1 — X) with (the free, U = 0, Green's function 
Go and) a self-energy X = E[/[G] given by a univer- 
sal functional. Apart from G, the universal functional 
Tiu must depend on the interaction parameters U only. 
Furthermore, the functional must satisfy a vanishing-curl 
condition or, alternatively, must be derivable from some 
(universal) functional <S>u[G] as THu[G] = 6<f>u[G]/6G 
(the temperature T is introduced for convenience). In 
short, "^-derivable" approximations are conserving. 

^-derivable approximations have been shown^ to ex- 
hibit several further advantageous properties in addition. 



One of these concerns the question of thcrmodynamical 
consistency. There are different ways to determine the 
grand potential of the system from the Green's function 
which do not necessarily yield the same result when us- 
ing approximate quantities. On the one hand, O may be 
calculated by integration of expectation values, accessi- 
ble by G, with respect to certain model parameters. For 
example, Q may be calculated by integration of the av- 
erage particle number, as obtained from the trace of G, 
with respect to the chemical potential \x. On the other 
hand, VL can be obtained as = $ + TrlnG - Tr(SG) 
without integration. A ^-derivable approximation con- 
sistently gives the same result for fi in both ways. 

At zero temperature T = there is another non-trivial 
theorem which is satisfied by any ^-derivable approxima- 
tion, namely Luttinger's sum rule££ This states that the 
volume in reciprocal space that is enclosed by the Fermi 
surface is equal to the average particle number. The orig- 
inal proof of the sum rule by Luttinger and Ward^ is 
based on the existence of <!> in the exact theory and is 
straightforwardly transferred to the case of a ^-derivable 
approximation. This also implies that other Fermi-liquid 
properties, such as the linear trend of the specific heat at 
low T and Fermi-liquid expressions for the T = charge 
and the spin susceptibility are respected by a ^-derivable 
approximation. 

There is a perturbation expansion^ which gives the 
Luttinger- Ward functional $t/[G] in terms of closed 
skeleton diagrams (see Fig. [1]). As a manageable $- 
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FIG. 1: Diagrammatic representation of the Luttinger- Ward 
functional $tr[G]. Double lines stand for the interacting one- 
particle Green's function G, dashed lines represent the ver- 
tices U . 
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derivable approximation must specify a (universal) func- 
tional $(7[G] that can be evaluated in practice, one usu- 
ally considers truncations of the expansion and sums up 
a certain subclass of skeleton diagrams only. This, how- 
ever, means that the construction of conserving approx- 
imations is restricted to the weak-coupling limit. 

One purpose of the present paper is to show that it is 
possible to construct ^-derivable approximations for lat- 
tice models of correlated fermions with local interactions 
which are non-perturbative, i.e. do not employ trunca- 
tions of the skeleton-diagram expansion. The idea is to 
employ the self-energy- functional theory (SFT) The 
SFT constructs the Luttinger-Ward functional $t/[G], or 
its Legendre transform Fjj [S] , in an indirect way, namely 
by making contact with an exactly solvable reference sys- 
tem. Thereby, the exact functional dependence of Fu [S] 
becomes available on a certain subspace of self-energies 
which is spanned by the self-energies generated by the 
reference system. 

The obvious question is whether those non- 
perturbative ^-derivable approximations have the 
same properties as the weak-coupling ^-derivable ap- 
proximations suggested by Baym and Kadanoff. This 
requires the discussion of the following points: 

(i) Macroscopic conservations laws. For fermionic lat- 
tice models, conservation of energy, particle number and 
spin have to be considered. Besides the static thermody- 
namics, the SFT concept concentrates on the one-particle 
excitations. For the approximate one-particle Green's 
function, however, it is actually simple to prove that the 
above conservation laws are respected. A short discus- 
sion is given in Appendix [X] 

(ii) Thermodynamical consistency. This issue has al- 
ready been addressed in Ref. [ll|. It has been shown that 
the [i derivative of the (approximate) SFT grand poten- 
tial (including a minus sign) equals the average particle 
number (N) as obtained by the trace of the (approxi- 
mate) Green's function. The same holds for any one- 
particle quantity coupling linearly via a parameter to the 
Hamiltonian, e.g. for the average total spin (S) coupling 
via a field of strength B. 

(hi) Luttinger sum rule. This is the main point to 
be discussed in the present paper. There are different 
open questions: First, it is straightforward to prove that 
weak-coupling ^-derivable approximations respect the 
sum rule as one can directly take over the proof for the 
exact theory. For approximations constructed within the 
SFT, a different proof has to be given. Second, it turns 
out that a non-perturbative <I>-dcrivable approximation 
respects the sum rule if and only if the sum rule holds 
for the reference system that is used within the SFT. 
As the original and thereby the related reference system 
may be studied in the strong-coupling regime, this raises 
the question which reference system does respect the sum 
rule, i.e. which approximation is consistent with the sum 
rule. Third, it will be particularly interesting to study 
reference systems which generate dynamical impurity ap- 
proximations (DIA)&£ and variational cluster approxi- 



mations (VCA) ) 10 i 12 as these consist of a finite number of 
degrees of freedom. Does the Luttinger sum rule hold for 
finite systems? Do the DIA and the VCA respect the sum 
rule? What is the simplest approximation consistent with 
the sum rule? Note that finite reference systems consist- 
ing of a few sites only have been shown£ i 13 i 14 i 15 i 16 i 17 i 18 i 19 
to generate approximations which qualitatively capture 
the main physics correctly. Finally, it is important to un- 
derstand these issues in order to understand whether and 
how a violation of the sum rule is possible within cluster 
extension a 20 i 21 i 22 ' 23 of the dynamical mean-field theory 
(DMFT)^^^ Note that the SFT com p r i S es the 

DMFT and certain^ 9 - cluster extensions and that possi- 
ble violations of the sum rule in the two-dimensional lat- 
tice models have been report ed ) 30 ' 31 i 32 including a study 
using the dynamical cluster approximation (DCA)^ 3 . 

The paper is organized as follows: A brief general dis- 
cussion of the Luttinger sum rule is given in the next 
section, and a form of the sum rule specific to systems 
with a finite number of spatial degrees of freedom is 
derived. Sec. [Till clarifies the status of the sum rule 
with respect to non-perturbative approximations gener- 
ated within the SFT framework. The results are elu- 
cidated by several numerical examples obtained for the 
most simple but non-trivial non-perturbative conserving 
approximation in Sec. lIVI Violations of the sum rule in fi- 
nite systems and their consequences are discussed in Sec. 
fVl Finally, Sec. I VII summarizes our main conclusions. 

II. LUTTINGER SUM RULE 

A system of interacting electrons on a lattice is gen- 
erally described by a Hamiltonian H(t,U) = Ho(t) + 
Hi(U) consisting of a one-particle part Ho(t) and an 
interaction Hi(U) with one-particle and interaction pa- 
rameters t and U, respectively. As a prototype, let 
us consider the single-band Hubbard mode l 34 i 35 ' 36 on a 
translationally invariant D dimensional lattice consist- 
ing of L sites with periodic boundary conditions. The 
Hamiltonian is given by: 

Here, i = 1, ...,L refers to the sites, a =1,1 is the spin 
projection, Ci„ (c| CT ) annihilates (creates) an electron in 
the one-electron state \io~), and n, CT = c^c^. Fourier 
transformation diagonalizes the hopping matrix t and 
yields the dispersion e(k). There are L allowed k points 
in the first Brillouin zone. 

Let G = Gt.u denote the one-electron Green's func- 
tion of the model H(t, U). In case of the Hubbard model, 
its elements are given by Gij(uj) = ((ci a ; c^))^. In the 
absence of spontaneous symmetry breaking, the Green's 
function is spin-independent and diagonal in reciprocal 
space. It can be written as Gk(uj) = l/(w + /i — e(k) — 
Sfc(w)) where y, is the chemical potential and Xfc(u;) the 
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sclf-cncrgy. Wc also introduce the notation St.u for 
the self-energy, and G t .o = 1/(uj + fJ, — t) for the free 
(non-interacting) Green's function which exhibits the de- 
pendence on the model parameters but suppresses the 
frequency dependence. Dyson's equation then reads as 
G t ,u = l/(G^o - £*,[/)■ 

The Luttinger sum rule££ states that 

(7V) = 2^6(G fc (0)) (2) 
k 



rule reads: 

2^^a m (fe)e(/i-o; m (fe)) 

k m 

= 2 E ( E - - E e ^ - ^»( fc )) J ■ 

k \ m n / 

(5) 

This form of the sum rule is convenient for the discussion 
of finite systems with L < oo. 



where N = J2ia nia * s the particle-number operator, 
(N) its (T = 0) expectation value, and the Heavy- 
side step function. The factor 2 accounts for the two 
spin directions. Since Gfc(O) -1 = fj, — e{k) — Sfc(0), the 
sum gives the number of k points enclosed by the in- 
teracting Fermi surface which, for L — > oo, is defined via 
fi— s(k) — Sfc(0) = 0. In the thermodynamic limit the sum 
rule therefore equates the average particle number with 
the Fermi-surface volume (apart from a factor (2tt) d /L). 
Note that, as 9(G fc (0)) = 9(1/G fe (0)), the sum rule Eq. 
@ also includes the so-called Luttinger volume^ which 
(for L — ► oo) is enclosed by the zeros of G&(0). 

The standard proof of the sum rule can be found in 
Ref. 1. It is based on diagrammatic perturbation theory 
to all orders which is used to construct the Luttinger- 
Ward functional $u [G] as the sum of renormalized closed 
skeleton diagrams (see Fig. [T]). We emphasize that the 
original proof straightforwardly extends also to finite sys- 
tems. For L < oo the sum in Eq. ^ is discrete. Actually, 
the proof is performed for finite L first, and the thermo- 
dynamic limit (if desired) can be taken in the end. The 
limit T — > 0, on the other hand, is essential and is re- 
sponsible for possible violations of the sum rule (see Sec. 

Below we need an alternative but equivalent formu- 
lation of the sum rule. We start from the following 
(Lehmann) representation for the Green's function: 

G»M=E ° m( *° m - (3) 

Here, ui m {k) — [i are the (real) poles and a m (k) the (real 
and positive) weights. For real frequencies uj, it is then 
easy to verify the identity: 

6(G fe (o;)) = Ee(w + M-w m (fe))-Ee(w + M-Cn(fe)) 

m n 

(4) 

where ( n (k) — A* is the n-th (real) zero of the Green's 
function, i.e. Gk{( n (k) — fx) = 0. 

For temperature T — we have (N) = 
2£fe f°^dw(-l/7r)ImG fe (w + i0+) and thus (N) = 
2 J2k Em a m(k)Q{n~uj m (k)). Hence, the Luttinger sum 



III. SELF-ENERGY-FUNCTIONAL THEORY 
AND LUTTINGER SUM RULE 

Within the self-energy- functional theory (SFT)^^^ 
the grand potential f2 is considered as a functional of the 
self-energy: 

n t , t7 [S]=Trln 1 - + (6) 

Here, the trace Tr of a quantity A is defined as TrA = 
T Hn 2 Efc e iu ^ )+ A k (iuj n ) where iuj n = i(2n + l)nT are 
the fermionic Matsubara frequencies, and the functional 
i*t/[X] is the Legendre transform of the Luttinger- Ward 
functional $[/[G]. The self-energy functional © is sta- 
tionary at the physical self-energy, JO^jy [E^jyJ/^X = 0, 
and, if evaluated at the physical self-energy, yields the 
physical value for the grand potential: £lt,u[Et.u] = 
Q tiU = -ThitT&q?{-0(H(t,U)- fiN)) where (3 = l/T. 
Comparing with the self-energy functional 

n t ,,u[n = Trln 1 - +F„\£\ (7) 

of a reference system with the same interaction but a 
modified one-particle part, i.e. with the Hamiltonian 
H(t', U), the not explicitly known but only [/-dependent 
functional Fjj [S] can be eliminated: 

f2 tit7 [£] = n t , it7 [E] + Trln 1 - Trln 1 . 

(8) 

An approximation is constructed by searching for a sta- 
tionary point of the self-energy functional on the sub- 
space of trial self-energies spanned by varying the one- 
particle parameters t': 

Inserting a trial sclf-cncrgy into Eq. © yields 

Hi, [/[£*', I/] = flt'.u + Trln— j — TrlnGf.Lr . 

(10) 

The decisive point is that the r.h.s. can be evaluated ex- 
actly for a reference system which is exactly solvable. 
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Apart from the free Green's function Gt.o, it involves 
quantities of the reference system only. 

This strategy to generate approximations has several 
advantages: (i) Contrary to the usual conserving approxi- 
mations, the exact functional form of Q^trfS] is retained. 
Any approximation is therefore non-pcrturbative by con- 
struction. On the level of one-particle excitations, macro- 
scopic conservation laws are respected as shown in Ap- 
pendixEJ (ii) With flt,u[^t' ,u] evaluated at the station- 
ary point t' = t' a , an approximate but explicit expression 
for a thcrmodynamical potential is provided. As all phys- 
ical quantities derive from this potential, the approxima- 
tion is thermodynamically consistent in itself (see Ref. 
[Til for details), (iii) As different reference systems gen- 
erate different approximations, the SFT provides a uni- 
fying framework that systematizes a class of "dynamic" 
approximations (see Refs. [29lji38l for a discussion). 

In the following we discuss the question whether or 
not a dynamic approximation respects the Luttinger sum 
rule. For this purpose consider hrst the Trln(- • • ) terms 
in Eq. (fTU|) . These can be evaluated using the ana- 
lytical and causal properties of the Green's functions 
as described in Ref. [9| (see Eq. (4) therein). Using 
-7Tn(l + exp(-w/T)) -> u9(-u) for T -> yields: 



Trln 



G t,0 ~ St'.U 



= 2^2^2(u) m (k) - n)Q{n - u> m (k)) 

fc m 

- 2 EZ)(Cn(fc)-A*)e(/i-Cn(fc)). 



(11) 



fc n 



Analogously, we have 



TrlnGW = 2£Xy m (JO-M)eO*-t4(*0) 

fc m 

- 2 EE(Cn(fc)-^)©(^-Cn(fc)). 

(12) 



fc n 



Note that the reference system is always assumed to be 
in the same macroscopic state as the original system, 
i.e. it is considered at the same temperature and, more 
importantly here, at the same chemical potential /x. Fur- 
thermore, it has been used that, by construction of the 
approximation, the self-energy and hence its poles at 
Cn(fc) — /i are the same for both, the original and the 
reference system. This implies that the second terms on 
the r.h.s. of Eq. (JTTJ) and iJHJ), respectively, cancel each 
other in Eq. (jTUjl . Finally, a (large but) finite system 
(L < oo) and a finite reference system are considered. 
Hence, the set of poles of the Green's function and of the 
self-energy as well as sums over k are discrete and finite. 

Taking the fx derivative on both sides of Eq. (fTU|) then 



yields: 

dSl t ,u[Vt>,u] 



an* 



t'U 



fc m 



2 X)E e ^ 



(13) 



fc rn 



Here we have assumed the ground state of the refer- 
ence system to be non-degenerate with respect to the 
particle number. From the (zero-temperature) Lehmann 
representation^ it is then obvious that, within a sub- 
space of fixed particle number, the /i-dependence of the 
Green's function is the same as its w-dependence, i.e. 
G(u>) — G(u> + fx) with a /i-indcpendent function G. Via 
the Dyson equation of the reference system, this prop- 
erty can also be inferred for the self-energy and, via the 
Dyson equation of the original system, for the (approx- 
imate) Green's function of the original system. Conse- 
quently, the poles of (G^q — Et^t/) -1 and of Gt>,u are 
linearly dependent on /i, i.e. uj m (k) and oj' m (k) in Eqs. 
(fTTj) and (fT2|) are independent of fj,. 

We once more exploit the fact that the self-energy of 
the original system is identified with the self-energy of the 
reference system. Using Eq. (U) one immediately arrives 
at 

(A) = (N)' + 2 J2 e(G fc (0)) - 2 J2 ©(GUO)) . (14) 
fc fc 

This is the final result: The Luttinger sum rule for 
the original system, Eq. @, is satisfied if and only if 
it is satisfied for the reference system, i.e. if (N)' = 
2£fc©(G fc (0)). 

A few remarks are in order. For the reference sys- 
tem, the status of the Luttinger sum rule is that of a 
general theorem (as long as the general proof is valid); 
(N)' and G' k (0) represent exact quantities. The above 
derivation shows that the theorem is "propagated" to the 
original system irrespective of the approximation that is 
constructed within the SFT. This propagation also works 
in the opposite direction. Namely, a possible violation 
of the exact sum rule for the reference system would im- 
ply a violation of the sum rule, expressed in terms of 
approximate quantities, for the original system. 

Eq. fTl")) holds for any choice of t'. Note, however, that 
stationarity with respect to the variational parameters i! 
is essential for the thermodynamical consistency of the 
approximation. In particular, consistency means that the 
average particle number (N) = — dflt,u[^t'.u]/dfj, on the 
l.h.s. can be obtained as the trace of the Green's function. 
Stationarity is thus necessary to get the sum rule in the 
form ([5]). 

There are no problems to take the thermodynamic 
limit (if desired) on both sides of Eq. (fl"4"]) (after divi- 
sion of both sides by the number of sites L) . The k sums 
turn into integrals over the unit cell of the reciprocal lat- 
tice. For a Z?-dimcnsional lattice the D — 1-dimcnsional 
manifolds of k points with Gfc(O) = oo or Gfc(O) = form 
Fermi or Luttinger surfaces, respectively. 
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For the above derivation, translational symmetry has 
been assumed for both, the original as well as the refer- 
ence system. Nothing, however, prevents us from repeat- 
ing the derivation in case of systems with reduced (or 
completely absent) translational symmetries. One sim- 
ply has to re-interprete the wave vector k as an index 
which, combined with to, refers to the elements of the 
diagonalized Green's function matrix G. The exact sum 
rule, Eq. (O, generalizes accordingly. The result (fl4)) re- 
mains valid (with the correct interpretation of k) for an 
original system with reduced translational symmetries. 
It is also valid for the case of a translationally symmet- 
ric original Hamiltonian where, due to the choice of a 
reference system with reduced translational symmetries, 
the symmetries of the (approximate) Green's function of 
the original system are (artificially) reduced. A typical 
example is the variational cluster approximation (VCA) 
where the reference system consists of isolated clusters of 
finite size. 



IV. TWO-SITE DYNAMICAL-IMPURITY 
APPROXIMATION 

While the Hartrec-Fock approximation may be con- 
sidered as the most simple weak-coupling ^-derivable 
approximation, the most simple non-pcrturbative $- 
derivable approximation is given by the dynamical- 
impurity approximation (DIA). This shall be demon- 
strated in the following for the single-band Hubbard 
model ([I]) as the original system to be investigated. The 
DIA is generated by a reference system consisting of a 
decoupled set of single-impurity Anderson models with 
a finite number of sites n s and is known^ to recover 
the dynamical mean-field theory in the limit n s — > oo. 
As long as the Luttinger sum rule holds for the single- 
impurity reference system, the DIA must yield a one- 
particle Green's function and a self-energy respecting the 
sum rule. 

The Hamiltonian of the reference system is H(t' , U) = 

Eiiffi with 

(7 a 
k=2 a k=2 a 

(15) 

For a homogeneous phase, the variational parameters 
t' = ({sq , , V^* }) can be assumed to be indepen- 
dent of the site index i: eq = £q, Ek = e£ , Vk = V^ l \ 
For the sake of simplicity, we consider the two-site DIA 
(n s = 2), i.e. a single bath site per correlated site only. 
In this case there are three independent variational pa- 
rameters only: the on-site energies of the correlated and 
of the bath site, Eq and e c = £k=2- respectively, as well as 
the hybridization strength V = Vk=2- As the reference 




FIG. 2: Filling dependence of the variational parameters 
at their respective optimized values and of the chemical po- 
tential. Calculations for the Hubbard model with a semi- 
elliptical free density of states of band width W — 4 and 
interaction strength U = W = 4 using the two-site DIA. 

system consists of replicated identical impurity models 
which are spatially decoupled, the trial self-energy is lo- 
cal and site-independent, Ejj(w) = SijT,(cu). 

Calculations have been performed for the Hub- 
bard model with a one-particle dispersion s(k) = 
Ylij e ~ ih ( Rt ~ Rs 'tij such that the density of one- 
particle energies D(e) is semi-elliptic. For |e| < W/2, 

D ( £ ) = jJ2 5 ( £ - £ ( fc )) = ^72 vW2) 2 - £ 2 • (16) 

k 

The free band width is set to W = 4. This serves as the 
energy scale. 

The computation of the SFT grand potential is per- 
formed as described in Ref. d. Stationary points of the 
resulting function Q(e ,£ c ,V) = £h,u\Ee ,e c ,v] are ob- 
tained via iterated linearizations of its gradient. There 
is a unique non-trivial stationary point (with F / 0). 
Fig. [U shows the variational parameters at this point as 
functions of the filling n. For the entire range of fillings, 
the ground state of the reference system lies in the in- 
variant subspace with N tot = J2a( c la c i<? + a l<j a i<?) = 2 - 
The parameters as well as the chemical potential are 
smooth functions of n. We have checked that the ther- 
modynamical consistency condition n = — L~ l dQ, / d = 

/-oo p( UJ )d LU is satisfied within numerical accuracy. Here 

p(w)=D(w + A*-S(w)) (17) 

is the interacting local density of states (DOS). 

At half-filling the values of the optimized on-site en- 
ergies are consistent with particle-hole symmetry. With 
Eo — (jl = —U/2 and e c — \i — the reference system is 
in the Kondo regime with a well-formed local moment at 
the correlated site. The finite hybridization strength V 
leads, for U = W, to a finite DOS p(ui = 0) > and thus 
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FIG. 3: Quasi-particle weight 2 as a function of the filling 
within the two-site DIA (full lines) and the two-site DMFT— 
(dashed lines). Calculations for U — W and U = 2W. 

to a metallic Fermi liquid as it is expected for the Hub- 
bard model within a (dynamical) mean-field description. 
Due to the simple structure of the self-energy generated 
by the two-site reference system, however, quasi-particlc 
damping effects arc missing. 

Decreasing the filling from n = 1 to n = drives the 
reference system more and more out of the Kondo regime. 
While e c stays close to the chemical potential, the on-site 
energy of the correlated site e crosses \x close to quar- 
ter filling and lies above \x eventually. Note that £0 = 
within the DMFT, i.e. for n s — > 00, while for finite n s 
there is a clear deviation from Sq = which is necessary 
to ensure thcrmodynamical consistency. For fillings very 
close to n = 0, the grand potential ftt,u[^e ,e c ,v] be- 
comes almost independent of S. This implies that it be- 
comes increasingly difficult to locate the stationary point 
with the numerical algorithm used. The slight upturn of 
£0 below n = 0.01 (see Fig. [5]) might be a numerical ar- 
tifact. 

It is instructive to compare the parameters with those 
of the two-site DMFT (2S-DMFT)42 The 2S-DMFT is 
a simplified version of the DMFT where a mapping onto 
the two-site single impurity Anderson model is achieved 
by means of a simplified self-consistency equation. As- 
suming £0 = as in the full DMFT, there are two pa- 
rameters left (e c and V) which are fixed by considering 
the first non-trivial order in the low- and in the high- 
frequency expansion of the self-energy and the Green's 
function in the DMFT self-consistency equation. Al- 
though being well motivated, this approximation is es- 
sentially ad hoc. One therefore has to expect that the 2S- 
DMFT is thermodynamically inconsistent and exhibits a 
violation of Luttinger's sum rule. A comparison of the 
DIA for n s = 2 with the 2S-DMFT is thus ideally suited 
to demonstrate the advantages gained by constructing 
approximations within the variational framework of the 
SFT. 

First of all, there are differences in fact. At half-filling 



the 2S-DMFT predicts the hybridization to be somewhat 
larger than the two-site DIA while the value for £ c is 
again fixed by particle-hole symmetry. Deviations grow 
with decreasing filling. Contrary to the two-site DIA, V 
monotonously increases and is larger in the entire filling 
range, £0 = by construction, and e c even diverges for 
n -> within the 2S-DMFT (see Ref. EJ. On the other 
hand, the system is essentially uncorrelated in the limit 
n — > 0. Strong differences in the parameters, which en- 
ter the self-energy only, therefore do not necessarily im- 
ply strongly different physical quantities. This is demon- 
strated by Fig. [3] which shows the quasi-particle weight 
calculated via 

as a function of the filling. While there are obvious differ- 
ences when comparing the results from the two-site DIA 
with those of the 2S-DMFT, the qualitative trend of z is 
very similar in both approximations. Both approxima- 
tions also compare well with the full DMFT: There is a 
quadratic behavior of z(n) for n — > 1 in the Fermi- liquid 
phase (U = W) and a linear trend when approaching the 
Mott phase (U = 2W). The critical interaction strength 
for the Mott transition is found to be U c « 1.46W for 
the two-site DIA and U c = 1.5W within the 2S-DMFT. 
For details on the Mott transition see Refs. Infill 

In case of a local and site-independent self-energy, the 
Luttinger sum rule can be written in the form^i 

H = no + E(w = 0) , (19) 

where /xq is the chemical potential of the free (U = 0) 
system at the same particle density. Eq. (fT9|) implies that 
not only the enclosed volume but also the shape of the 
Fermi surface remains unchanged when switching on the 
interaction. Using Eq. (|17[) this immediately implies^ 

p(0) = D(jm>) = po(0) , (20) 

i.e., in case of a correlated metal, the value of the inter- 
acting local density of states at u) = is independent of 
U and thus fixed to the value of the density of states of 
the non-interacting system at the same filling. 

The interacting and the non-interacting DOS arc plot- 
ted in Fig. 2] for different fillings and for U = W and 
U = 2W. The impurity self-energy of the two-site ref- 
erence system is an analytical function of ui except for 
two first-order poles on the real axis. Via Eq. (|17p this 
two-pole structure implies that the DOS consists of three 
peaks the form of which is essentially given by the non- 
interacting DOS. At half-filling the three peaks are eas- 
ily identified as the lower and the upper Hubbard band 
and the quasi-particle resonance as it is characteristic 
for a (dynamical) mean-field description.— For U = W 
the resonance still has a significant weight. The weight 
decreases upon approaching the critical interaction, and 
the resonance has disappeared in the Mott insulator for 
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FIG. 4: Interacting local density of states p{uS) (solid lines) 
for different fillings as indicated. Calculations using the two- 
site DIA for U = W (left) and U = 2W (right). For n = 0.25, 
n — 0.5 and n = 0.75 the non-interacting DOS pa{u) is shown 
for comparison (dashed lines). Note that p(0) = /9o(0). The 
dotted line for U = 2W in the top panel is the DOS for 
n = 0.99. 



U = 2W. Hole doping of the Mott insulator is accom- 
plished by the reappearance of the resonance at u) = 
which preempts the creation of holes in the lower Hub- 
bard band42 As can be seen in the spectrum for n = 0.99 
in the top panel (dotted line), the quasi-particle res- 
onance appears within the Mott-Hubbard gap. With 
decreasing filling, the upper Hubbard band gradually 
shifts to higher excitation energies and loses weight. This 
weight is transfered to the low-energy part of the spec- 
trum. For lower fillings where the Kondo regime has been 
left, one would actually expect that the quasi-particle res- 
onance disappears by merging with the lower Hubbard 
band. This, however, cannot be described with the sim- 
ple two-pole structure of the self-energy. One therefore 
should intcrprcte the gap around = — 1 at n = 0.25 
as an artifact of the approximation. Furthermore, the 
widths of the Hubbard bands is considerably underes- 
timated as damping effects are missing completely. The 
filling-dependent spectral-weight transfer across the Hub- 
bard gap as well as the energy positions of the main 
peaks, however, are in overall agreement with general 
expectations! 34 ' 43 

It is worth emphasizing that this simple two-site 
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FIG. 5: Numerical results for the difference between the 
volume enclosed by the Fermi surface Vfs and the filling n as 
a function of n for U = W = 4. The Luttinger sum rule (Vfs — 
n = 0) is exactly respected by the two-site DIA. Results for 
the 2S-DMFT and the Hubbard-I approximation are shown 
for comparison. Dashed line: difference between the filling n 
and the average occupation of the correlated (impurity) site 
in the reference system at stationarity for the two-site DIA. 



dynamical-impurity approximation exactly fulfills the 
Luttinger sum rule. In Fig.[¥]this can be seen by compar- 
ing with the DOS of the non-interacting system (dashed 
lines). The non- interacting DOS cuts the interacting one 
at u> = which shows that Eq. (f20|) is satisfied. Note 
that this is trivial for n = 1 as this is already enforced 
by particle-hole symmetry. Off half-filling, however, the 
pinning of the DOS to its non-interacting value at u = 
is a consequence of $-derivability and thereby a highly 
non-trivial feature. 

In contrast, the 2S-DMFT docs show a violation of 
Luttinger's sum rule which, however, must be attributed 
to the ad hoc nature of the approximation. Fig. [5] shows 




0.2 0.4 0.6 0.8 1 
filling n 

FIG. 6: Filling dependence of the compressibility k for U = 
W as obtained within the 2S-DMFT via k = dn/dfi (solid 
line) and via a general Fermi- liquid relation (Eq. \22\ . dashed 
line). Using the two-site DIA identical results are obtained 
for both cases. 
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the difference between the volume enclosed by the Fermi 
surface 

2 r° 
V FS = - V 9(/i-£(fe)-E(0)) = 2 / deD(e+n-Z(0)) 

L k J ~™ 

(21) 

and the filling n as a function of the filling. As can be 
seen, there is an artificial violation of the sum rule for the 
2S-DMFT which is of the order of a few per cent while 
for the ^-derivable two-site DIA the sum rule is fully re- 
spected. Note that, unlike the DMFT and also unlike 
the simplified 2S-DMFT, the two-site DIA predicts a fill- 
ing which slightly differs from the average occupation of 
the correlated impurity site in the reference system (see 
dashed line in Fig. [5]) . For a finite number of bath sites n s 
this appears to be necessary to fulfill the Luttinger sum 
rule. The figure also shows the result obtained within 
the Hubbard-I approximation^ Here a very strong (ar- 
tificial) violation of up to 100 % (for n close to half-filling) 
is obtained. This should be considered as a strong draw- 
back which is typical for uncontrolled mean-field approx- 
imations. 

There are more relations which, analogously to the 
Luttinger sum rule, can be derived by means of per- 
turbation theory to all orders^ in the exact theory and 
which are respected by weak-coupling conserving approx- 
imations. For example, the compressibility, defined as 
K = dn/dfj,, can be shown to be related to the interact- 
ing DOS and the self-energy at the Fermi edge via 

K = M0) (22, 

Fig. [5] shows that for the 2S-DMFT it makes a difference 
whether n is calculated as the /^-derivative of the filling 
or via Eq. ([22]) . Again, this must be attributed to the 
fact that the 2S-DMFT is not a <I>-derivable approxima- 
tion. Contrary, the two-site DIA does respect the general 
Fermi-liquid property (|22p and thus yields the same re- 
sult in both cases (see Fig.[6j). 

V. VIOLATION OF LUTTINGER'S SUM RULE 
IN FINITE SYSTEMS 

The preceding section has demonstrated that the two- 
site DIA satisfies the Luttinger sum rule. According to 
Eq. (fT4|) , we can conclude that the Luttinger sum rule 
must hold for the corresponding reference system, i.e. for 
the two-site single-impurity Anderson model. Of course, 
this can be verified more directly by evaluating Eq. (O. 
In case of a finite system or a system with reduced trans- 
lational symmetries, the Green's function is a matrix with 
elements G a p{uj) where a refers to the one-particle basis 
states, and the Luttinger sum rule reads: 

]T a«6( M - U W) = J2 ^ "ft) - E - ■ 

k,m k,m k,n 

(23) 



Here the index k labels the elements of the diagonalized 
Green's function, i.e. Eq. ([5]) is generalized by replacing 
(fc, a) — > k. In case of an impurity model, Eq. (f23|) ac- 
tually represents the Friedel sum rulei 44 ' 45 For the two- 
site single-impurity Anderson model, the different one- 
particle excitation energies — /i, the zeros of the 
Green's function — M and the weights a^ 1 are eas- 
ily determined by full diagonalization. We find that Eq. 
(f2"3"| is satisfied in the entire parameter space (except for 
V = 0, see below). 

Note that a violation of the sum rule occurs when, as 
a function of a model parameter x, a zero of the Green's 
function crosses u> = for x = x c . At x c the number of 
negative zeros counted by the second term on the r.h.s. 
changes by one while the first term as well as the l.h.s. 
remain constant since (unlike a pole) a zero of the Green's 
function is generically not connected with a change of the 
ground state (level crossing). This implies that the sum 
rule would be violated for x < x c or for x > x c . 

The case V = is exceptional. Within the two-site 
DIA this corresponds to the Mott insulator (see Fig. 2J 
topmost panel for U = 2W). For V = the reference 
system consists of two decoupled sites, and the Green's 
function becomes diagonal in the site index. There is 
no zero of the local Green's function corresponding to 
the uncorrelated site. We can thus concentrate on the 
correlated site where the local Green's function exhibits 
a zero at r] — /i = Eq + U/2. In the sector with one electron 
at the correlated site (so < [i < £o + U), the second 
term on the r.h.s. changes by two at /i = /i c = Eq + U/2 
because of the two-fold degenerate ground state. In this 
case Luttinger's sum rule in the form (j2"3"|) is violated for 
/i < fx c and for fi > fi c . This "violation", however, is a 
trivial one which immediately disappears if the ground- 
state degeneracy is lifted by applying a weak field term, 
for example. 

Fig. [7] shows a phase diagram of the single-impurity 
Anderson model with n s = 4 sites as obtained by 
full diagonalization. The diagram covers the entire 
range of the total particle number N = J2a( c t c v) + 
Ha E r kU( a L a k*) from N = to N = 2n s = 8. A non- 
degenerate ground state is enforced by applying a small 
but finite magnetic field. No violation of the Luttinger 
sum rule is found. We have repeated the same calculation 
also for n s = 10 using the Lanczos technique*^ Again, 
the sum rule is found to be always satisfied (We have 
performed calculations for different U and bath parame- 
ters). This might have been expected as the (n s — > oo) 
Anderson model can generally be classified as a (local) 
Fermi liquid4£ 

The situation is less clear in the case of correlated lat- 
tice models such as the Hubbard or the t-J model. For 
two dimensions there are several numerical studies using 
high-temperature expansion^ quantum Monte-Carlo^i 
extended DMFT , 32 ' 48 and dynamical cluster approxima- 
tion (DCA)^ 3 - which indicate a violation in the strongly 
correlated metallic phase close to half-filling. For studies 
of large clusters or studies directly working in the thcrmo- 
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FIG. 7: Phase diagram /i vs. e of the single-impurity An- 
derson model with n s = 4 sites. Total particle numbers are 
indicated by Roman figures. Results have been obtained by 
full diagonalization for the following model parameters. One- 
particle energies: eq = (correlated site), eu = £ + (k — 3) 
with k — 2,3,4 (uncorrelated bath sites). Hubbard interac- 
tion: U = 2e. Hybridization strength: Vk = 0.1 for k — 2, 3, 4. 
To lift Kramers degeneracy in case of an odd particle number, 
a weak (ferromagnetic) field of strength b = 0.001 is coupled 
to the local spins. The dashed line marks the particle-hole 
symmetric case. Luttinger's sum rule is found to be satisfied 
in the entire parameter space. 



dynamic limit, a definite conclusion on the validity of the 
sum rule is difficult to obtain as finite-temperature or ar- 
tificial broadening effects etc. must be controlled numer- 
ically. Contrary, full diagonalization of Hubbard clusters 
consisting of a few sites only can provide exact results. 
While their direct relevance for the thermodynamic limit 
is less clear, it is important to note that reference sys- 
tems with a finite number of sites or a finite number of 
correlated sites provide the basis for a number of cluster 
approaches within the SFT framework. Via Eq. (| 14[) their 
properties are transferred to the approximate treatment 
of lattice models in the thermodynamic limit. 

The validity of Eq. (|23[) has been checked for Hubbard 
clusters of different size and in different geometries. The 
(i vs. U phase diagram for an L = 4-sitc open Hubbard 
chain with nearest-neighbor hopping in Fig. [H] shows a 
representative example. Again, a small but finite field 
term is added to avoid a ground-state degeneracy. As the 
chemical potential, for fixed U, is moved off the particle- 
hole symmetric point fi — U/2 and exceeds certain criti- 
cal values (red lines), the particle number N [as obtained 
from the l.h.s. of Eq. (|23|) ] changes from N = L down to 
(up to) N = (N = 2L). A critical /i value indicates a 
change of the ground state (level crossing) that is accom- 
panied by a change of the ground-state particle number. 
In the one-particle Green's function this is characterized 

(k) 

by a pole wV — [i crossing lo — 0. The blue lines indicate 
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FIG. 8: Phase diagram fi vs. U of the Hubbard model with 
L = 4 sites (open chain) as obtained by full diagonalization. 
Nearest-neighbor hopping t = — 1. A weak (ferromagnetic) 
field of strength 6 = 0.01 is applied to lift Kramers degener- 
acy. The dashed line marks the particle-hole symmetric case. 
Particle numbers (l.h.s. of Eq. ([23])) are indicated by Roman 
figures. R.h.s. of Eq. (J23J : Arabic figures. Luttinger's sum 
rule is found to be violated for sufficiently strong U. U c i, U c .2- 
critical interactions. 



those chemical potentials at which a zero of the Green's 

(k) 

function Q — fi crosses u) = 0. Whenever this happens 
the r.h.s. of Eq. (|23|) changes while the l.h.s. is constant. 
Fig. [5] shows that this occurs several times in the N = L 
sector. At the particle-hole symmetric point \i = U/2 
the Luttinger sum rule is obeyed while it is violated in a 
wide region of the parameter space corresponding to half- 
filling N = L. However, a critical interaction strength 
U c turns out to be necessary. The value for U c strongly 
varies for different cluster sizes and geometries but has 
always been found to be positive and finite. Note that 
for L = 4 the sum rule is fulfilled for any particle num- 
ber N ^ L. Qualitatively similar results can be found for 
the L = 2-sitc Hubbard cluster where calculations can be 
done even analytically. Again, a violation of the sum rule 
is found in the half-filled sector beyond a certain critical 
U. 

This has already been noticed by Rosch^ and was used 
in combination with a strong-coupling expansion to ar- 
gue that a violation of the sum rule generically occurs 
for a Mott insulator. Stancscu et aliSl have shown quite 
generally that the sum rule is fulfilled when particle-hole 
symmetry is present (the Luttinger surface is the same 
as the Fermi surface of the non-interacting system) but 
violated in the Mott insulator away from particle-hole 
symmetry. It is interesting to note that these arguments 
cannot be used to construct a violation of the sum rule 
within DMFT or for a single-impurity Anderson model: 
For an (almost) particle-hole symmetric case and model 
parameters describing a Mott insulator (within DMFT), 
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FIG. 9: Ground-state particle number (red Roman figures, 
l.h.s. of Eq. (I23[l ) and prediction by the Luttinger sum rule 
(blue Arabic figures, r.h.s. of Eq. (J23J) as functions of the 
chemical potential for a L = 9-site Hubbard cluster with pe- 
riodic boundary conditions. Arabic numbers are only given 
when different from Roman ones. Calculations using the 
Lanczos method and a finite but small magnetic field and 
finite but small on-site potentials to lift ground-state degen- 
eracies. 



an odd number of sites n s (with n s — > oc) must be consid- 
ered and thus a magnetic field is needed to lift Kramers 
degeneracy. Even an infinitesimal field, however, leads 
(at zero temperature) to a finite and even large polar- 
ization corresponding a well-formed but unscreened local 
moment. This polarization is incomplete for any finite 
U as the DMFT predicts a small but finite double occu- 
pancy for a Mott insulator. Still there is a proximity to 
the fully polarized band insulator which finally results in 
a weakly correlated state and thus in a situation which 
is unlikely to show a violation of the sum rule. 

We have also considered Hubbard clusters with L = 9 
and L = 10 sites by using the Lanczos technique^ 
Calculations have been performed for different Lanczos 
depths Z max to ensure that the results are independent 
of Zmax- Fig- M displays an example for L = 9 and a 
highly symmetric cluster geometry with periodic bound- 
ary conditions and a well-defined reciprocal space. To lift 
ground-state degeneracies resulting from spatial symme- 
tries as well as the Kramers degeneracy, small but fi- 
nite on-site potentials and a small magnetic-field term 
are included in the cluster Hamiltonian. Fig. [TO] shows 
an example for L = 10 sites without any spatial sym- 
metries. Kramers degeneracy for odd N is removed by 
applying a small magnetic field. With the figures we com- 
pare the expressions on the left-hand and the right-hand 
side of Eq. (j2"5)) . Obviously, the sum rule is respected in 
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FIG. 10: The same as Fig. H but for 10 sites. 



most cases. Violations are seen for half-filling N = L, 
i.e. in the "Mott-insulating phase", which is consistent 
with Ref. H^. However, the sum rule is also violated 
in the "metallic phase" close to half-filling, namely for 
N = L — 1 (Fig. [9J L = 9) and N = L - 1, L - 2 (Fig. 
flOl L = 10). This nicely corresponds to the generally 
observed tren d 30 ' 31 i 32 ' 33 i 48 for violations in the slightly 
doped metallic regime. We have also verified that the 
sum rule is restored by lowering U. 

Fig . l9l and [TOl demonstrate that the sum rule is violated 
in the whole /i range corresponding to TV = L — 1. This 
is an important point as it shows that it is irrelevant 
whether the T = limit is approached by holding (N) 
fixed and adjusting /i = fi(T) or by fixing \i and let (N) = 
(N)(T) be T-dependcnt. A violation of the sum rule is 
found in both cases. 

Kokalj and Prelovsek^i have demonstrated that viola- 
tions of the sum rule can also be found for the t-J model 
on a finite number of sites. Our result provides an ex- 
plicit example showing that not only for t-J— but also 
for Hubbard clusters a violation can be found when the 
chemical potential is set to /i = limj^o fi(T) with (a(T) 
obtained for given (N) = const. Anyway, the original 
proofs does not depend on this choice for /i but appears 
to work for any \x. 

The results raise the question which assumptions used 
in the original proof of the theorem arc violated or where 
the proof breaks down. Note that the recently proposed 
alternative topological proof^ assumes a Fermi-liquid 
state from the very beginning and thus cannot be applied 
to a finite system. Using weak symmetry-breaking fields, 
a more or less trivial breakdown due to ground-state de- 
generacy has been excluded. An analysis of the ground 
state of the L = 2 and L = 4 Hubbard clusters which 
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are accessible with exact (analytical or numerical) meth- 
ods has shown that, for model parameters where the sum 
rule is violated, the interacting ground state can never- 
theless be adiabatically connected to the non-interacting 
one. This excludes level crossing as a potential cause for 
the breakdown. While we cannot make a definite state- 
ment, it appears at least plausible that the violation of 
the sum rule results from a non-commutativity of two 
limiting processes, the infinite skeleton-diagram expan- 
sion and the limit T — > 0. 

Using a functional-integral formalism, the Luttinger- 
Ward functional at finite T can also be constructed in a 
non-perturbative way, i.e. avoiding an infinite summation 
of diagrams, as has been shown recently— Formally, the 
Luttinger sum rule can be obtained by exploiting a gauge 
invariance of the Luttinger- Ward functional [see Ref. [53| : 

-l—* u [G(iw n ))=0. (24) 
o(iu> n ) 

If at all, this invariance can only be shown for T = 
where iu> n becomes a continuous variable. Unfortunately, 
the non-perturbative construction of $>jj requires a T > 
formalism. Hence, the validity of the sum rule depends 
on question whether the limit T — > commutes with the 
frequency differentiation. Necessary and sufficient condi- 
tions for this assumption are not easily worked out. An 
understanding of the main reason for the possible break- 
down of the sum rule in finite systems, very similar to 
the case of Mott insulators, is therefore not yet available 
(see also the discussion in Ref. l49h . 

VI. CONCLUSIONS 

^-derivable approximations are conserving, thermody- 
namically consistent and, for T = 0, formally respect cer- 
tain non-trivial theorems such as the Luttinger sum rule. 
As the construction of the Luttinger- Ward functional <& 
is by no means trivial and may conflict with the limit 
T — ¥ or different other limiting processes, however, the 
validity of the sum rule may be questioned. Violations of 
the sum rule can be found in fact for the case of strongly 
correlated electron systems. For Mott insulators and fi- 
nite systems in particular, a breakdown is documented 
easily. 

This implies that a general approximation for the 
spectrum of one-particle excitations (of the one-particle 
Green's function) may violate the sum rule for two pos- 
sible reasons, namely because (i) the sum rule is violated 
in the exact theory, or (ii) the approximation generates 
an artificial violation. 

Within the usual weak-coupling conserving approxima- 
tions, such as the fluctuation-exchange approximation, 
the sum rule always holds as the formal steps in the gen- 
eral proof of the sum rule can be carried over to the 
approximation - but with the important simplification 
of a limited class of diagrams. This also implies that 
weak-coupling conserving approximations, when applied 



beyond the weak-coupling regime, might erroneously pre- 
dict the sum rule to hold. 

The present paper has focussed on non-perturbative 
conserving approximations. Non-perturbative approxi- 
mations, constructed within the framework of the self- 
energy-functional theory and referring to a certain ref- 
erence system, are ^-derivable and consequently respect 
certain macroscopic conservation laws and are thermo- 
dynamically consistent. Whether or not the sum rule 
holds within the approximate approach, however, cannot 
be answered generally. We found that Luttinger's sum 
rule holds within an (SFT) approximation if and only if 
it holds exactly in the corresponding reference system. 

The reference system that leads to the most simple 
but non-trivial example for a non-perturbative conserv- 
ing approximation consists of a single correlated and 
a single bath site. For this two-site system, we have 
found the sum rule to be valid in the entire parameter 
space. Consequently, the resulting two-site dynamical- 
impurity approximation (DIA) - opposed to more ad hoc 
approaches like the two-site DMFT - fully respects the 
sum rule as could be demonstrated in different ways. In 
view of the simplicity of the approximation this is a re- 
markable result. Since the sum rule dictates the low- 
frequency behavior of the one-particle Green's function, 
important mean-field concepts, such as the emergence of 
a quasi-particle resonance at the Fermi edge, are quali- 
tatively captured correctly, even away from the particle- 
hole symmetric case. This qualifies the two-site DIA for a 
quick but rough estimate of mean-field physics, including 
phases with spontaneously broken symmetries. 

Full diagonalization and the Lanczos method have 
been employed to show that also the single-impurity An- 
derson model with a finite number of n s > 2 sites respects 
the sum rule. Consequently, this property is transferred 
to an n s -site DIA. For n s — ► oo the full dynamical mean- 
field theory is recovered which is thereby recognized as 
the prototypical non-perturbative conserving approxima- 
tion. Clearly, in the case of the DMFT, $-derivability is 
well known^i and obvious, for example, when construct- 
ing the DMFT with the help of the skeleton-diagram ex- 
pansion. 

Using as a trial self-energy the self-energy of a cluster 
with L > 1 correlated sites, generates an approximation 
where short-range spatial correlations are included up to 
the cluster extension. These variational cluster approxi- 
mations provide a first step beyond the mean-field con- 
cept. Again, whether or not the sum rule is respected 
within the VCA depends on the reference system itself. 
For the L — 2 Hubbard cluster, analytical calculations 
straightforwardly show that violations of the sum rule 
occur at half-filling, beyond a certain critical interaction 
strength. In the thermodynamic limit, this would corre- 
spond to the Mott-insulating regime. Applying the Lanc- 
zos method to larger clusters, has shown, however, that 
a breakdown of the Luttinger sum rule is also possible 
for fillings off half-filling. For sufficiently strong U, the 
sum rule is violated in the whole N = L — 1-particlc 
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sector. This would correspond to a (strongly correlated) 
metallic state in the thermodynamic limit. Whether or 
not a VCA calculation is consistent with the sum rule, 
then depends on the set of cluster hopping parameters t' 
which make the self-energy functional stationary. First 
VCA calculations 5 -^ for the D = 2 Hubbard model at low 
doping and using clusters with up to L = 10 sites do 
predict a violation in fact. 

It is by no means clear a priori what happens in a clus- 
ter approach using additional bath degrees of freedom as 
variational parameters, as e.g. in the cellular DMFT . 21 ' 22 
The usual periodization of the self-consistent C-DMFT 
self-energy, however, should be avoided when testing the 
sum rule as this introduces an additional (though physi- 
cally motivated) approximation. Instead. Eq. {5) must be 
used with k re-interpreted as an index referring to the el- 
ements of the self-consistent diagonalizcd lattice Green's 
function. 

Employing the dynamical cluster approximation 
(DCA)2£ represents an alternative which directly oper- 
ates in reciprocal space. From a real-space perspective, 
the DCA is equivalent with the cellular DMFT but ap- 
plied to a modified model H = H(t, U) — > H (t, U) with 
modified hopping parameters which are invariant under 
supcrlatticc translations as well as under translations on 
the cluster! 29 ! 55 In the limit L — > oo the replacement 
t — ► t becomes irrelevant. Analogous to the C-DMFT, 
the sum rule then holds within the DCA if and only if 
it holds for the individual cluster at self-consistently de- 
termined cluster parameters. Note, however, that this 
requires that (besides the DCA self-energy) the modified 
hopping t instead of the physical hopping has to be con- 
sidered in the computation of the volume enclosed by 
the Fermi (Luttingcr) surface of the lattice model. This 
is exactly what is usually done in DCA calculations. 

Within this context and in view of the violations found 
for finite Hubbard clusters, it is possible to understand 
why a non-perturbative cluster approximation, like the 
VCA,^i or a cluster extension of the DMFT, like the 
DCA^ can produce results that are inconsistent with 
Luttinger's theorem. 
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APPENDIX A: MACROSCOPIC CONSERVATION 
OF ENERGY, PARTICLE NUMBER AND SPIN 

The one-particle Green's function as obtained within 
an approximation generated by the choice of a reference 
system respects the macroscopic conservation laws which 



result from symmetries of the system with respect to con- 
tinuous transformation groups: 

Energy conservation is apparently respected as by con- 
struction the approximate SFT Green's function depends 
on a single frequency only, i.e. is invariant under time 
translations. 

Conservation of the total particle number and spin is 
respected if the approximate G transforms in the same 
way as the exact Green's function under global U(l) and 
SU(2) gauge transformations. Consider a general trans- 
formation of the form 

4-ct=E*v«4 ( A1 ) 



with unitary S such that the interaction part Hi (U) 
of the Hamiltonian is invariant (a refers to the states 
of the one-particle basis). In a diagrammatic approach, 
the invariance of H\(U) implies that the corresponding 
conservation law is respected "locally" at each vertex. 
Hence, for a conserving approximation in the sense of 
Baym and Kadanoff, the transformation behavior of the 
free Green's function is then propagated by the diagram 
rules to the full Green's function. Consequently, the lat- 
ter must transform under S in the same way as the exact 
G, i.e. 

G af) ->&*,, = (SGSfl) afl . (A2) 

Consider now the case of the SFT. One has to show 
that the approximate Green's function G for the trans- 
formed system with Hamiltonian H is given by G — 
SGS^ if G is the approximate Green's function of the 
model H. Applying the transformation (|A1|) to H, one 
finds H = H (t) + Hi (U) —>Tl= H (t) + Hi(U) with 
t = StS^ . Again, S is assumed to leave the interaction 
part invariant. 

The Green's function G of the transformed model is 
(approximately) constructed via 



from the free Green's function of the transformed model 
and the SFT self-energy which is the self-energy of the 
reference system H' = Ho(t s ) + Hi(U) at the stationary 
point t s . 

For the transformed problem H, the stationary point 
t s is determined from the SFT Euler equation: 




(A4) 

As an ansatz to solve the Euler equation we take 

t = St'iS^ (A5) 

with t[ to be determined. The transformation law (|A2[) 
for the exact Green's function of the reference system 
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is Gji v = Gg t i ,u = SGt^.u ■ This also holds for 
the free Green's function. Using the Dyson equation of 
the reference system we can deduce v = SUf yrS^. 
Furthermore, for the free Green's function of the trans- 
formed original model we have Gj = SGt.oS^ . Using 
these results, we see that Eq. (|A4[) is equivalent to 

^[g^-^u' Gt '- U ) dt> -°- 

(A6) 

But this is just the Euler equation for the original 



model which is solved by t[ = t' s . Remembering the 
ansatz made, we now have for the stationary point 
t' s = St? a £fl. Inserting this into Eq. ([A3]) gives G = 
S(G^l - S t -, t/ )- 1 S' t = SGS* which is the desired re- 
sult. 

Note that one has to ensure that the stationary point 
for the transformed problem t s = St'gS'' lies within the 
space of one-particle parameters characteristic for the ref- 
erence system. For models with local interaction part 
and for local (and also global) gauge transformations, 
however, this is always easily satisfied. 
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